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Abstract 

We consider several problems that involve finding the eigenvalues and 
generating the eigenstates of unknown unitary gates. We first examine 
ControUed-U gates that act on qubits and assume that we know the eigen- 
values. It is then shown how to use singlet states to produce qubits in the 
eigenstates of the gate. We then remove the assumption that we know 
the eigenvalues and show how to both find the eigenvalues and produce 
qubits in the eigenstates. Finally, we look at the case where the unitary 
operation acts on qutrits and it has two eigenvalues of 1 and one of —1. 
We are able to use a singlet state to produce a qutrit in the eigenstate 
corresponding to the —1 eigenvalue. 



1 Introduction 

A common problem, which arises in quantum mechanics, is finding the eigen- 
values and eigenstates of an operator, usually the Hamiltonian. The eigenvalues 
are the values that the observable corresponding to the operator can assume, 
and the eigenstates are the states of the system in which that observable will 
have a definite value. 
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With the advent of quantum algorithms, a natural question to raise is 
whether there are quantum algorithms that will efficiently find the eigenval- 
ues and eigenvectors of operators. The answer to this question is, in fact, yes. 
Based on earlier work by Kitaev Cleve, et al. developed an algorithm that 
can estimate an eigenvalue if one copy of the eigenstate is provided This 
algorithm was analyzed further by Abrams and Lloyd They pointed out 
that it is not necessary to have a copy of an eigenstate to use this procedure. 
One can start with an arbitrary input state, and at the end of the procedure, 
one will obtain an eigenvalue corresponding to an eigenstate that has a nonzero 
overlap with the input state. This is not a deterministic procedure; we could 
obtain any eigenvalue whose eigenstate has a significant overlap with the input 
vector. In addition, Abrams and Lloyd showed that at the output one has not 
only the eigenvalue, but a set of qubits that is in a state that is a good approx- 
imation to the eigenstate corresponding to the measured eigenvalue. How good 
this approximation is was recently investigated by Travaglione and Milburn Q . 

This procedure requires that one have some knowledge about the eigenstate 
one is trying to generate and whose eigenvalue one is trying to find. In particular, 
it is necessary that the input state have a substantial overlap with the desired 
eigenstate. It may be possible to accomplish this if there is some information 
available that allows a guess for the state to be made. For example, when 
finding the ground state energy of a not-too-complicated Hamiltonian, it might 
be possible on physical grounds to obtain a rough idea of what the ground state 
would look like, and this information could be used to design an appropriate 
input state for the eigenvalue estimation algorithm. In many cases, however, 
there will be little if any information to guide ones choice, with the result that 
the input state may have a very small or no overlap with the desired eigenstate. 

Here we shall show that a different input state, a singlet state, allows one to 
find all eigenvalues and eigenstates of an unknown ControUed-U gate simulta- 
neously. We shall start with the case of a single qubit gate where we know the 
eigenvalues and wish to generate output qubits in the eigenstates of the gate. 
We shall then proceed to the case where the gate still operates on only a single 
qubit, but we do not know either its eigenvalues or its eigenstates. Our object 
then is to find the eigenvalues and produce output qubits in the eigenstates. 
Finally, we shall consider a ControUed-U gate that acts on qutrits and has two 
eigenvalues of 1 and one eigenvalue of —1. It is possible to use a singlet state 
to produce a qutrit in the eigenstate corresponding to the eigenvalue —1 with 
a network that contains only two ControUed-U gates. This procedure is easily 
generalized to D-dimensional quantum systems, qudits. Given a ControUed-U 
gate that acts on qudits and has eigenvalues 1, which is {D — l)-fold degenerate, 
and —1, it is possible to produce a qudit in the eigenstate with eigenvalue — 1 
by using a network containing D — 1 ControUcd-U gates. 
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2 Generation of eigenstates of a Controlled- ?7 
gate with known eigenvalues 

Consider the following problem. We are given a Controlled-?/ gate which acts 
on single qubits and wc would like to generate its eigenstates. We know that 
the eigenvalues of the gate are 1 and —1, but have no information about its 
eigenstates. A measurement-based strategy for doing this would involve sending 
qubits through the gate and measuring them. For example, we could use the 
basis states |0) and |1) to obtain information about the matrix elements of U in 
that basis. If we send the state |l)a|0)b, where a is the control bit and b is the 
target bit, through the gate, the probability that target qubit at the output is 
in the state |0)& can be measured. This probability is just equal to |(0|C/|0)p. 
The probability that the target qubit is in the state is just |(l|?7|0)p. These 
measurements give us information about two of the matrix elements of U, and 
the fact that we know that the eigenvalues are 1 and —1 means that these are 
the only two we have to know. In particular we have that 

(0|C/|0) = -(l|f/|l) (0|f/|l) = (l|f/|0)*. (1) 

Information about the relative phase of these matrix elements can be gained 
by using the input state (|0) + e*^|l))/-\/2. The probability, po that the output 
vector is in the state |0) is 

Po = mU\0)+e^O{0\U\l)\'. (2) 

After sending through many qubits we would have an estimate of the matrix 
elements, and we could then diagonalize the matrix. This information could 
then be used to generate qubits in the eigenstates of U. This procedure involves 
many qubits and many uses of the gate. What we shall now present is a quantum 
strategy that will produce both eigenstates with certainty using only 3 qubits 
and requiring only one use of the gate. 

Let a and b be the control and target qubits of the gate, as before. To this 
we add a third qubit, which we shall denote by c. We now define the following 
states 

\±x) = ^(I0)±|l)) 

|<A,) = -L(|01)-|10)), (3) 

and note that the rotational invariance of the singlet state, \4>s), implies that it 
can also be expressed as 

\cl>,) = ^{\u+)\u.)-\u.)\u+)), (4) 

where \u^) is the eigenstate of U with eigenvalue 1 and is the eigenstate 
with eigenvalue —1. More explicitly, the transformation specified by 

y|o) = |«+) v\i) = |«_), (5) 
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is unitary, and, because \(j>s) is invariant under U{2) (X) (7(2), we have 



{V^VMs)^\^s). (6) 

Because the phase of the eigenstates is arbitrary, in general Eqs. (Q) and (||) wiU 
be true only up to an overall phase factor, but we shall assume that the phases 
of and |u_) are such that the phase factor is equal to one. 
We now start our 3-qubit system in the state 

|*m)afcc = I + X)a\(f>s)bc- (7) 

After qubits a and b go through the ControUed-C/ gate the state of the system 
is 

\'^out)abc = ^[\0)a{\u+)b\U-)c ~ \U-)b\u+)c) 

+ \l)ai\u+)b\u^)c + \U-)b\u+)c)]. (8) 

We now want to measure the a qubit in the | ± x) basis. In order to see the 
result of doing so we can express |5'o«t)afcc as 

\'^out)abc = -^(1 + X)a\u+)b\u^)c - \ - x) a | W-) f, | U+) c) • (9) 

This equation implies that if we measure a and get | + x) , then qubit b is in the 
1 eigenstate and qubit c is in the —1 eigenstate, while if we get \ — x), then it 
is just the other way around. Therefore, with one use of the gate we have, with 
certainty, generated its two eigenstates. 

Now let us see what happens when the eigenvalues are not ±1 but two other 
complex numbers of modulus unity. We shall denote the eigenstates as 
and \u2) and the corresponding eigenvalues as as e^^^ and e*^^ respectively. We 
again use the same 3-qubit scheme and choose the input state to be \+x)a\(t>s)bc- 
The output state is now given by 

\'^out)abc = —q{\vi)a\ui)b\u2)c - \v2) a\u2) b\ui) c) , (10) 

where 

\vi) = -^(|0)+e^^Ml» 

h) = -i=(|0)+e^«^|l)). (11) 

Note that \vi) and \v2) are not orthogonal, but because we assume we know 
the eigenvalues, these vectors are known. At this point we can apply the op- 
timal procedure for distinguishing two nonorthogonal states [^-[0. This is a 
generalized measurement that can be applied when we are given a state, which 
is one of two known states, and we want to determine which of the two states 
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it is. The measurement has one of three possible outcomes; it cither tchs us 
without error which of the two states we have, or it teUs us that it has failed 
to distinguish the states. Applied to and \v2), the procedure would succeed 
with a probability of 

1 - \{vi\v2)\ = 1 - -^[1 +cos(ei - 6*2)]'/^ (12) 

The closer the phase difference between the eigenvalues is to tt, the greater the 
probability of success of this procedure. If \v-i) is detected at output a, then 
\ui) is at output b and |u2) at c. On the other hand, if \v2) is detected at a, 
then \u2) is at b and is at c. 

If one allows more than one use of the gate there are other possibilities. 
Suppose that we know the eigenvalues of U are 1 and i. Then the eigenvalues 
of U"^ are 1 and —1. We can then apply the above procedure to generate the 
eigenstates with one modification. The three-qubit initial state is the same, but 
qubits a and b pass through two Controlled-/7 gates instead of one. Again qubit 
a is measured at the output in the | ± x) basis. If | + a;) is found, then the 
eigenstate corresponding to 1 is at output b and that corresponding to i is at 
output c. If I — a;) is found, the b and c outputs are reversed. 

As preparation for the next section, let us consider one last, harder problem. 
Suppose that we know that the eigenvalues are not the same, and that they 
are members of the set {1,-1, i, —i}. In this case we have partial rather than 
complete information about the eigenvalues, and we again want to generate 
qubits in the eigenstates. This can be done using four qubits, a Controlled-?/ 
gate, and a Controlled- gate. This last gate can be constructed from two 
Controlled-[/ gates in sequence. Qubit a is the control bit for the Controlled- 

gate, qubit b is the control bit for the Controlled- [/ gate, and qubit c is the 
target bit for both. Let and [112) be the eigenstates of U with corresponding 
eigenvalues Zi and Z2, where zi and Z2 are members of the set {1,-1, 
The input state is 

|*m)afecd = \+x)a\+x)b\(j>s)cd 

= \+x)a\+x)b^-^=j {\ui)c\u2)d-\U2)c\ui)d)- (13) 

The output state of the network is given by 

\^out)abcd = Gac{U'^)Gbc{U)\^in)abcd 

= ^[(|00)ab + Zi\01)ab + zf\W)ab + zf\U) ab)\ui) c\u2) d 

-(|00)„6 + Z2|01)afc + zl\10)ab + z| 1 1 1) afc) I W2)c ^1 )d] , (W) 

where Gjk{U^) is the operator corresponding to a Controlled- C/" gate with 
control bit j and target bit k. Now consider the vector 

\v{z))ab = limab + Z\0l}ab + Z^\lO)ab + ^'|ll)a6). (15) 
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The vectors \r]{l))ab, |»7(-l))ah, \v{i))ab, and \r]{-i))ab form an orthonormal 
basis of the space of the two qubits a and b. If we now measure the two qubit 
system ab in this basis, we can determine one of the eigenvalues of U, e.g. if 
the result of the measurement is |77(l))ab, then one of the eigenvalues is 1. The 
eigenstate corresponding to this eigenvalue will emerge from output c, and the 
eigenstate corresponding to the other, unknown, eigenvalue will emerge from 
output d. 

This procedure will allow us to find one of the eigenvalues, if we know that 
they belong to a hmited set, and generate both eigenstates. A better proce- 
dure would allow us to find both eigenvalues, remove the restriction that they 
belong to a particular set, and generate both eigenstates. Such an algorithm is 
presented in the next section. 

3 Application of phase estimation to find un- 
known eigenvalues and eigenvectors 

Suppose that we have an unkown ControUed-U gate, and we want to find its 
eigenvalues and generate qubits in its eigenstates. This can be done by modi- 
fying the phase estimation algorithm of Cleve, et al. and using a singlet state 
as the input [||. One takes two phase-estimation circuits for the same gate, 
and sends into each one of two particles, which, together, form a singlet state. 
This avoids the main disadvantage of the original algorithm. There, besides the 
ControUed-U gates, one also needed a qubit prepared in one of the eigenstates. 
Sending this qubit through the network would then generate an estimate of the 
eigenvalue for this eigenstate. An alternative is to send in a random qubit, in 
which case one gets an estimate for a random eigenvalue. In particular, the es- 
timate corresponds to one of the eigenvalues whose eigenstates have a nonzero 
overlap with the input state. The original qubit is left in a state that is close 
approximation to the eigenstate corresponding to the measured eigenvalue Q . 

Each of the two networks, which we shall label A and B is constructed as fol- 
lows. We have n control qubits, which for network A we shall call Al, A2, . . . An, 
and one target bit, which we shall call A. Each of the control bits is initially in 
the state ( 1 0) -I- 1 1 ) ) / a/2. Control bit Aj is connected to a gate that does nothing 
if the control bit is 0, and performs the operation U^^ if the control bit is 1. 
The network B is identical. The effect of the entire network is given by 

|*out) = G(B„)b(C/'"") . . .G(bi)b({/)G(^„)a([/2"") . . . G(^l)A(C/)|*m). 

(16) 

Let the eigenstates of U be \ui) and \u2) with eigenvalues e*"^^ and e*"^^, respec- 
tively. As before, the singlet state can be expressed in terms of these eigenstates 

\(I)s)ab = -^(|ui)a|w2)b - \u2)a\ui)b)- (17) 
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The initial state of the system is then 



n-l 



n-1 



= -1-{\u,}a\u2}b - \U2}A\m)B)llmAj + \l}Aj) UHO) Bk + \l}Bk). 

(18) 



j=0 fe=0 

After passing through the networks this state becomes 
1 



|*,n) 



v/22" 



n-1 



n-1 



j=0 



k=0 



n-1 



n-1 



\U2)Al[(\^)Aj+e^''^^\l)Aj)\Ui)Bl[(\^)Bk + e'' *N1>B.) 
j=0 fc=0 



.(19) 



The products in the above equation can be expressed as a sum over n-digit 
binary numbers. For example, 



2"-l 



j=0 y=0 



(20) 



The first digit of the n-digit binary number y corresponds to the state of system 
An, the second to that of A{n — 1), and so on. In the above equation we have 
indicated this cxpUcitly with subscribts on the; state, but if the future these 
will be omitted and this correspondence will be understood. It is still necessary 
to indiciate whether \y) is a state of An, . . .Al or Bn, . . .Bl, and this will be 
indicated by the subscripts A and B, respectively. We then have that 



I ^initial) 



\/22" 



2"-l 2"-l 

\ui)a{Y1 e'^^y\yh)\u2)B{Y,e^^n^h) 

y—0 w=0 



(21) 



The next step is to apply the quantum inverse Fourier transform operation 
to states A and B. This takes the state \y) to 



(22) 



Before applying this, however, we want to express the phases 4>i and 02 in 
different way. First, let xj = (pj /{2n) for j = 1, 2, which implies that < Xj < 1. 
In addition, let Xj be the closest integer to (we assume Xj is expressed in 
binary form) so that 

^i = ^+<5,-, (23) 
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where \Sj\ < 1/2"+^. If we now apply the inverse Fourier transform wc find that 

2"/2" 



2"-l 



g2^iy[(3?i/2")+<5j]|^^ _^ ^ ^ ^ ^ g27rjy(x,-z)/2"g27riy5j|^^ 



„27j-iy[(3;,/2")+<5j] 

i/=o ~ v=a z=o 

It is possible to perform the y summation in the above equation 

I I — g27ri5j2" 

^ 2" 1 - e2"[((^j-^)/2")+'53] ■ 



(24) 



(25) 



This function is peaked about z = Xj and the maximum value of its magnitude 
is greater than 2/n [Q. After applying the inverse Fourier transform to both 
systems A and B, our state is 



1 



2"-l 2"-l 

*i>^(X! 5(2;^i,<5i)|z)^)|M2)i3( X! 9is;X2,S2)\s)g) 

2"-l 2"-l 



\u2)a{Y^ g{z;x2,S2)\z)^)\ui)BiY^ gis;xi,Si)\s)g) 

z=0 s=0 



(26) 



We now measure both systems A and B in the computational basis. The 
most likely results are either Xi for A and X2 for B, in which case qubit A is 
in I Ml) and B is in |u2), or X2 for A and Si for B, in which case qubit A is in 
|w2) and B is in In either case we have both the eigenvalues (to n places 

in base 2) and qubits in the eigenvectors. 



4 Higher dimensional systems 

The reasoning in the previous sections can be extended from qubits to qu- 
dits, D-dimensional quantum systems. The fully antisymmetric state of D D- 
dimensional quantum systems is a J7 (D) singlet . If we denote the computa- 
tional basis states by \n), where n = 0, 1, ... 13 — 1, this state can be expressed 
as 

D-l D-1 

\MD)) = 7= E ■ • ■ E ^n-M • • • \Jd), (27) 

where is the totally antisymmetric tensor of rank D. Now consider a 

unitary operator U whose eigenstates are \uj), where j = 1, . . . D. The fact that 
\(j>s{D)) is a singlet means that it can be expressed as 

D 

^^i-JuWji) ■■■\ujn), (28) 

jn = l 



\MD)) = 
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where e*'' is a phase factor that depends on how the phases of the oigcnstates 
are chosen. We shah subsequently assume that they have been chosen so that 
M = 0. 

Let us now consider the following problem for the case Z) = 3; its general- 
ization to the case of arbitrary dimension is straightforward. We are given a 
Controlled-J7 gate, where the control is a qubit and the target is a qutrit. If 
the control qubit is in the state |0) nothing happens to the target qutrit, and 
if it is in the state |1), the operation U is performed on the qutrit. This gate 
corresponds to the operator Gjk{U), where j is the control qubit and k is the 
target qutrit. We shall assume that the operator U has eigenvalues 1 and — 1, 
where the eigenvalue 1 is degenerate, and we would like to produce a qutrit in 
the eigenstate corresponding to —1. 

This can be done with two Controlled-f7 gates, two qubits, and three qutrits. 
The initial state of the system is 

|*m)a...e = | + x)c,\ + x)b\M^))cde- (29) 

Particles a and b arc qubits and c, rf, and e arc qutrits. Qubit a is the control 
bit for qutrit c and qubit b is the control bit for qutrit d. The output state is 
given by 

|*o«t)a...e = Gac{U)GUU)\^in)a...e- (30) 

Let \ui) and \u2) be orthonormal cigcnstatcs of U with eigenvalue 1 and \v) be 
the eigenstate with eigenvalue —1. In terms of these states we have that 

\^out)a...e = "^[1 ~ + x) b{\vUiU2) cde - \vU2Ui)cde) 

+ 1 + X)a\ - x)b{\u2VUi)cde - \uiVU2)cde) 

+ 1 + X)a\ + x)b{\uiU2V)cde " \U2U1V) cde)]- (31) 

We now measure qubits a and b in the | ± x) basis. If we find a in the | — x) 
state, then qutrit c is in the eigenstate with eigenvalue —1, while if qubit b is 
in the | — x) state, then it is qutrit d that is in the —1 eigenstate. Finally, if 
both of these qubits is found to be in the | + x) state, then qutrit e is in the —1 

eigenstate. 

In the corresponding problem for qudits, U has eigenvalues —1, which is 
nondegenerate, and 1, which is D — 1 fold degenerate. The object is to produce 
a qudit in the eigenstate corresponding to —1. To do so one uses a network 
consisting oi D — 1 qubits, D qudits in a singlet state, and D — 1 Controlled-J7 
gates. The procedure is a simple generalization of the one just discussed for 
qutrits. 

5 Conclusion 

We have shown that singlet states in combination with Controlled- J7 gates can 
be used to produce qubit, or qudits, in eigenstates of the operator U. If U is the 
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evolution operator corresponding to some Hamiltonian, its eigenstates are just 
those of the Hamiltonian. This procedure will not tell us what those eigenstates 
are, but we can perform measurements on the qudits in those states in order to 
gain information about them. We may also simply be interested in performing 
further operations on these states without measuring them first, and we now 
have a way of producing them. 

If singlet states are combined with the phase-estimation algorithm for finding 
eigenvalues of J7, we can, in a certain sense, diagonalize the operator. We saw 
that for qubits we both knew the eigenvalues, at least to a level of approximation 
that we can determine, and we produced qubits in states that are very close to 
the eigenstates of U. This procedure should generalize to qudits. 
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